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The Kardar-Parisi-Zhang (KPZ) equation with infinitesimal surface tension, dynamically develops 
sharply connected valley structures within which the height derivative is not continuous. We discuss 
the intermittency issue in the problem of stationary state forced KPZ equation in 1 + 1-dimensions. 
It is proved that the moments of height increments C a =< \h(xi) — h{x2)\ a > behave as \xi — X2^ a 
with £ a = a for length scales \xi — X2\ « o. The length scale a is the characteristic length of the 
forcing term. We have checked the analytical results by direct numerical simulation. 
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I. INTRODUCTION 

The growth, formation and morphology of interfaces 
has been one of the recent interesting fields of study be- 
cause of its high technical and rich theoretical advan- 
tages. On account of the disorder nature embedded in 
the surface growth, stochastic differential equations have 
been used as a suitable tool for understanding the behav- 
ior of various growth processes. Such equations typically 
describe the interfaces at large length scales, which means 
that the short length scale details has been neglected in 
order to derive a continuum equation by focusing on the 
coarse grained properties. A great deal of recent theo- 
retical modeling has been started with the work of Ed- 
ward and Wilkinson describing the dynamics of height 
fluctuations by a simple linear stochastic equation [1-5]. 
By adding a new term proportional to the square of the 
height gradient, Kardar, Parisi and Zhang made an ap- 
propriate description for lateral interface growth [6] . The 
(1 + l)-dimcnsional forced KPZ equation is written as 



h t {x, t) = —(h x ) 2 + uh xx + f(x, t) 



(1) 



where a > and f(x,t) is a zero-mean, statistically ho- 
mogeneous, white in time random force with covariance 



(f(x, t)f(x',t')) = 2D D(x - x')S(t - t'). 



(2) 



Typically the spatial correlation of the forcing is con- 
sidered to be a delta function, mimicking the short length 
correlation. Here the spatial correlation is considered as 



D(x-x') = 



1 



(trj 2 ) 1 ^ 



exp( ) 



(3) 



where a is much less than the system size L, i.e. a « L, 
which represents a short range character for the random 
forcing. The KPZ equation has come to famous the 
"Ising model" of non-equilibrium physics. It is indeed 
the simplest equation nevertheless capturing the main 



determinants of the growth dynamics i.e. nonlincarity, 
stochasticity, and locality. The theoretical richness of 
the KPZ model is partly due to close relationships with 
other areas of statistical physics. It is shown that there 
is a mapping between the equilibrium statistical mechan- 
ics of a two dimensional smcctic-A liquid crystal onto 
the non-equilibrium dynamics of the (1+1)- dimensional 
stochastic KPZ equation [7]. It has been shown in [8] 
that, one can map the kinetics of the annihilation pro- 
cess A + B — > with driven diffusion onto the (1+1)- 
dimensional KPZ equation. Also the KPZ equation is 
closely related to the dynamics of a sine-Gordon chain 
[9], the driven-diffusion equation [10,11], high T c - super- 
conductor [12] and directed paths in the random media 
[13-26] and charge density waves [27], dislocations in dis- 
ordered solids [3], formation of large-scale structure in 
the universe [28-31] , Burgers turbulence [32-60] and etc. 

It is useful to rescale the KPZ equation as h! = h/h , 
r' = r/r and t' = t/t . If we let h = (-^) 1/2 and 

r 2 

to = where rg is a characteristic length, all of the 
parameters can be eliminated, except the coupling con- 
stant g = a Ji ■ The limit g — > oo (or zero tension limit, 
v — >• 0), is known as the strong coupling limit [60]. Phase 
diagram information extracted from the renormalisation 
group flow indicates that d = 2 plays the role of a lower 
critical dimension. For d < 2, the Gaussian fixed point 
(a = 0) is infrared-unstable, and there is a crossover 
to the stable strong coupling fixed point. For d > 2, a 
third fixed point exists, which represents the roughening 
transition. It is unstable and appears between the Gaus- 
sian and strong coupling fixed points which are now both 
stable. Only the critical indices of the strong-coupling 
regime (g — > oo or v — > 0) are known in 1+1 dimen- 
sions and their values in higher dimensions as well as 
properties of the roughening transition have been known 
only numerically [61-67], and the various approximation 
schemes [68-76]. 

For finite a, in the strong coupling limit {y — » 0) non- 



linear term in the KPZ equation will dominate. The non- 
linearity of the KPZ equation in this limit includes the 
possibility of singularity formations in a finite time as a 
result of the local minima instability. Meaning that there 
is a competition between the diffusion smoothing effect 
( the Laplacian term), and the enhancement of non-zero 
slopes. In one spatial dimension the sharp valleys are de- 
veloped in a finite time. As indicated in figs.(l) and (2), 
the geometrical picture consists of a collection of sharp 
valleys intervening a scries of hills in the stationary state 
[77]. 

The main difficulty with the KPZ equation is that it 
is controlled, in all dimensions, by a strong disorder ( or 
strong coupling) fixed point and efficient tools are miss- 
ing to calculate the exponents and other universal prop- 
erties e.g. scaling functions, amplitudes, etc. Despite the 
fact that in one dimension, the exponents are known, but 
many properties, including the probability density func- 
tion (PDF) of the height of a growing interface have been 
so far measured only in numerical simulations. 

In this article, the statistical properties of the KPZ 
equation in the strong coupling limit [y — > 0) is investi- 
gated. The limit is singular, i.e. through which the sur- 
face develops sharp valleys. So starting with a flat surface 
after a finite time scale, t c ~ (7r) 1 / 6 J D ( r 1 a" 2/3 er 5/3 [77], 
sharp valley singularities are dynamically developed. In 
the singular points spatial derivative of the field h{x, t) 
is not continuous. One of the main problem in this area 
is the scaling behavior of moments of height increments 
C a =< \h(xi) — h{x2)\ a > and the probability density 
function (PDF) of 5h = h(xi) - h(x 2 ), that is P(Sh). 
Inspired by the methods proposed recently by Weinan E 
and Vanden Eijnden [47], a statistical method is devel- 
oped to describe the moments of the height and height 
gradient increments. We derive a master equation for 
joint PDF of the height and its gradient increments in 
1+1 dimensions. It is shown that in the stationary state 
where the sharp valleys arc fully developed, the relax- 
ation term with infinitesimal surface tension leading to 
an unclosed term in the PDF's equation. However we 
show that the unclosed term can be expressed in terms 
of statistics of some quantities defined on the sharp val- 
leys. We identify each sharp valley in position ?/o with 
three quantities, namely the gradients of h in the posi- 
tions y + 1 2/o - an d hs height from the h. The dynamics 
of these quantities are given in [77]. Here it is proved 
that to leading order, when \x± — a^l << a i fluctua- 
tion of the height field is not intermittent. The ana- 
lytic form of the amplitudes of the structure functions is 
also given. The absence of the intermittency means that 
C a =< \h(xi) — h(x2)\ a > scales as \xi — x 2 \^ a , where 
£ a is a linear function of a. It is proved that for length 
scales \x\ — x 2 \ « <J, the exponents £ a are equal to a. 

The paper is organized as follows; in section two, 
the known results for the moments of height increments 
C a =< \h(xi) — h(x 2 )\ a > ,for length scales \xi — x 2 \ » 



a, are expressed. In section three, we derive the master 
equation for the joint PDF of height and its gradient in- 
crements for given surface tension v and for length scales 
\x\ — x 2 \ « a. It is shown that the surface tension term 
makes the PDF's equation unclosed. In section 4 we will 
consider the limit of v — > of the master equation and 
derive the scaling exponents of height increments mo- 
ments. Also a comparison between the analytical results 
and direct numerical simulation are given. Details of cal- 
culations are presented in the appendices A and B. 



II. SCALING EXPONENTS OF 
HEIGHT-DIFFERENCE MOMENTS FOR 
FORCED KPZ EQUATION AND FOR THE 
LENGTH SCALES |Xi - X 2 \ » a. 

In this section a review of the known results for the 
scaling exponent of height increments moments for the 
KPZ equation in 1 + 1-dimcnsions with white in time 
and space forcing is given. Indeed the limit a — > is 
considered in equation (2). In this limit the equation (2) 
can be written as follows; 



{f{x, t)f{x'.t')) = 2D 5(x - x')5(t - f ) 



(4) 



For this type of forcing, U[h(x, t)], the probability func- 
tional of h(x,t) = h(x,t) — (h) satisfies the functional 
Fockcr-Planck equation [1,2], 



d d x 



-J—[(%(Vh) 2 +uV 2 h)ll} 
6h{x) K 2 ' 



+D / d d x 



-n, 



Sh 2 {x) 

where its solution in the (l+l)-dimensions is 



n = cxp [ 



2,Dn 



dx(h x y 



(5) 



(6) 



Therefore if one introduce G{x — x') = (h(x)h(x')) as 
a Green's function, then it satisfies the following differ- 
ential equation 



d xx G{x — x') = -S(x — x') 



(7) 



so that (h(x)h(0)) = —^f-\x\. Now we can write the 
second moment of height increments for small it's as 



(\h(x) - h(0)\ 2 ) = —\x\ 



(8) 



In a similar way it can be seen that the higher mo- 
ments, (\h(x) — h(0)\ a ) scale with x as \x\f , which means 
that for the length scale a << \x\ — x 2 \ « L, the expo- 
nents arc £ a = |. 



There are a few comments on the result obtained for 
the functional PDF, equation (6). It is evident that the 
probability density functional (in 1 + 1-dimensions ) is 
independent of the coefficient of the nonlinear term i.e. 
a, so the result is independent of the strength of the cou- 
pling constant. As it can be seen the scaling relation 
is similar to an ordinary random walk problem. If one 
considers the random force with smooth spatial correla- 
tion, the problem changes to a more complicated one and 
there is no any closed solution for the functional PDF. In 
the next sections we will show that the moments of the 
height increments for the length scales \xi — x 2 \ << <J has 
the scaling exponents £ a = a and the amplitude of the 
moments are depend on the coefficient of the nonlinear 
term a . 



III. THE MASTER EQUATION GOVERNING 
THE PROBABILITY DENSITY FUNCTION OF 
THE HEIGHT-DIFFERENCE AND 
GRADIENT-DIFFERENCE FOR GIVEN 
SURFACE TENSION 

In this section , focusing on the (l+l)-dimcnsional 
KPZ equation and it's corresponding Burgers equation, 
the master equation describing the evolution of the joint 
two point PDF, P{h{x\) — h{x2),u(xi) — u{x%)) of the 
height and corresponding height gradients increments is 
derived. The (1+1) dimensional KPZ equation is written 



h t(x, t) = -h 2 x + vh xx + f(x, t) 



(9) 



where a > and f(x,t) is a zero-mean, statistically ho- 
mogeneous, white in time random force. Its covariancc 
is given by cq.(2). Using the map —d x h = u, the corre- 
sponding Burgers equation is written as 



ut = -auu x + vu xx - f x (x, t). 



(10) 



Defining the two point generating function as 
A2, fii, H2, Xi, X2, t) = (0), where is defined as, 



:= cxp (— iXhi — iXhz — ifJ-iUi — ^2^2)- 



(11) 



The fields h x and h 2 are the height of the surface 
at points x\ and x 2 . The fields u\ = —d Xl hi and 
1*2 = —d X2 h 2 are related to the corresponding height 
gradients. As it is seen the generating function is the 
ensemble average of 0. The time evolution of Z will be 



Z t = -iXi(h lt O) - i\ 2 (h 2t &) 
- i^i(u lt Q) - i/i 2 (w 2 t@}- 



(12) 



Using the equations (9) and (10) and noting that, in 
equation (12), hi, h 2 , u\ and U2 can be substituted by 



h } -» *aI7>> -» ^b' Ul l WI and U2 the 
time evolution of Z can be rewritten as 



7 - aAl /<"v\ , - a ^ 2 /r\\ 

d d 

- a/xi- — (ui Xl Q) - a[i 2 - — (u 2x2 Q) 

Ofil 0[X2 

- iAi(/ie) -a 2 (/ 2 e) + m(h Xl e) + in 2 (f 2xa e) 

+ i\iv(u lxi Q) + i\ 2 v(u 2x2 Q) - i(jbiv(ui xx Q) 

- i^ 2 v(u 2xx Q). (13) 
Now using , 

(u jx Q) = — (©>*,. + — 3 = 1,2 (14) 

Mi /'.; 

the equation governing Z can be written as 

Z t = l — (©Wl + l — ("Wa 

-^^(1(6)^+^(9^) 

V\ll ^1 

-a^ 2 — ( — (0).2+A 2 (0) Al2 ) 

OfJ.2 ^2 

+ <e> Xl +A 1 (0) Ml ) 
Mi 

+ i\ 2 v{— (0),, +A 2 (0) M2 )+^ + £. (15) 

M2 

Here J- and Q stand for 

J r =-iAi(/ie)-*Aa(/ 2 e> 



+ iiJ>i(fixi@) + im{f-2x 2 Q) 
Q = ifi :>' " 1 , , + i^2v(u 2xx Q). 



(16) 



In equation (15) the terms J 7 and Q are the only terms 
which arc not closed respect to Z. Indeed the term T 
can be also closed according to Novikov's theorem , 

T = ( - (A? + X 2 2 )K{0) - 2AiA 2 A»)Z 

+ (-(/"? + (4)K(0) - 2 t M W2 K(x)) Z (17) 

where K(x) = 2DqD(x) and x = Xi — x 2 . So G is the 
only term preventing equation (15) to be closed which 
can be referred to a sort of dissipative anomaly. 

The PDF P(hi,h2,ui,U2,xi,X2,t) is defined as the 
two-point joint probability density function (PDF) at the 
points Xi and X2 with their related heights hi and h 2 , and 
their gradients ui and u 2 . The PDF can be constructed 
by Fourier transforming the generating function Z 



, . f dXi dX 2 dpi dfj.2 

P(h 1> h 2 ,u 1> u 2 ,x 1) x 2 ,t) = J — — — — 

x exp (iXhi + iXh 2 + ifiiui + i\i 2 u 2 )Z{X\, X 2 , fn, n 2 , Xi,x 2 , 

Fourier transformation of equation (15) gives the fol- 
lowing equation for the PDF 

d d a d d d 2 

dui du 2 * 2 dh\ du\ du 2 1 

_*_d__d__d_ ( 2 P) 

2dh 2 d Ul du 2 y 2 ' 

- a-^—P x + a^—P x 
au 2 ou\ 

d d . 
ohi au 2 

Q Q 

dui du 2 ^ 1 x 

d d . 

+ a a — "5 — \U2Px) 
am ou 2 



(18) 



where x = x 2 — xi, y = Xl + X2 ; d Xl = —d x + \d y and 
d X2 = d x + \d y . The terms Z{[i\\i 2 T} and C{^\^ 2 Q) are 
the Fourier transformations of equations (16) and (17), 
multiplied in [ii and fx 2 , where for L{\i\\i 2 T) is 



d d , d 2 



i) 2 



2^, d d d d p 
du\ du 2 dhi dh 2 



8 d , d 2 



i:) 2 



k xx ®) q Ui q U2 (q u 2 du 2 ^ 

d 2 d 2 
duf OU2 



(19) 



and C(fiifi 2 G) is defined as 

CfaimiG) := -u{{ui Xlxl \hi,h 2 ,ui,u 2 ,x)P} u ^ u2 

- v{{u 2x2X2 \hi,h2,u 1 ,U2,x)P} uiU2U2 . (20) 
For later use we define G as 
G := G(hi, h 2 , ui,u 2 ,x,t) + G(hi, h 2 , u\, u 2 , -x,t) (21) 
where 

G(hi,Ui,Ui,U2,x,t) = -v{(ui XlXl \hi,h 2 ,ui,u 2 ,x)P} ui . 

Also we can simply substitute C(fiifi 2 Q) with 
G UlU2 . In cq.(20), (u lxiXl \h 1 ,h 2 ,u 1 ,u 2 ,x) and 



(u 2x2X2 \hi, h 2 , u\, 1x2, x) are the averages of U\ XlXl and 
U2x 2 x 2 conditional that the heights and velocities fields be 
tfii, I12, u\ and U2 with a spatial difference x. Now we are 
interested in writing an evolution equation for the PDF's 
of height and its gradients difference. We change the 
variables h\,h 2 , u\ and u 2 with u\ = u — u 2 = u + % , 
hi = h — I and h 2 = h + f . Integrating over u and h 
the PDF of the height and height gradient difference is 
obtained 

P S (H,u,x,t) 



I 



dhduP(h-i,h+^,u-^,u+^,x,t) (22) 



Finally using the eq.(18), the master equation can be 
written as, 



P^t = ~2aPL - a{wP 5 x ) uw + 2(fc(0) - k(x))P^ 



2(k xx (0)-k xx (x))PL^ + G: 



6 



(23) 



where by considering the definition of G in (21), G 6 
would be 



G 5 {^,u,x,t) 



dhduG. 



(24) 



It is clear that the G s , which is proportional to sur- 
face tension v, makes the master equation unclosed. In 
appendix A, it is proved that for finite a in the limit of 
v — > 0, G s can be written in terms of the quantities which 
arc defined on singularities. 



IV. CALCULATION OF THE MOMENTS AND 
NUMERICAL SIMULATION 

As shown in the previous section, the presence of sur- 
face tension is, makes the master equation unclosed. 
However in the limit v — > 0, that is the KPZ equa- 
tion with an infinitesimal surface tension, one can 
find exact scaling exponents of the moments (\h(xi) — 
h{x2)\ a \u(x\) — u(x 2 )\ b ). It should be noted that the 
M-field satisfying the Burgers equation, for finite cr's, de- 
velops discontinues or shock solutions in the limit v — ► 0. 
Consequently for finite a the height field shows up as a 
set of sharp valleys at the positions where the shocks are 
located, continuously connected by some hill configura- 
tions, as indicated in Figs.(l) and (2). As mentioned, 
each sharp valley in position yo is identified by three 
quantities, namely the gradients of h in positions y + , 
y - and it's height from h. It is evident that the term 
lim, y ^o vu xx is zero at the positions where no sharp valley 
exists. Therefore in the limit v — ► 0, only small intervals 
around the sharp valleys will contribute to the integral 
in the eq.(24). Within these intervals, boundary layer 
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FIG. 1. Different snapshots of the time-evolution of the 
height, with correlation lengths a ~ L/10 for the random 
periodic force, until the time that the system finally reaches 
to it's stationary state. The average distance between the 
sharp valleys is of order of a [56,57]. 




FIG. 2. Different snapshots of the time-evolution of the 
height, with correlation lengths a ~ L/100 for the random 
periodic force, until the time that the system finally reaches 
to it's stationary state. 
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FIG. 3. In the upper graph the log-log plot of 
< \h(0) — h(x)\ a > vs x, for moments a = 0.7,5 and 7 are 
demonstrated. In the middle figure the corresponding scaling 
exponent £ a for height increments are plotted. The £ a has a 
linear dependence on a. In the lower figure the scaling expo- 
nent £j, for the moments of the height gradients increments 
are shown. 



analysis can be used for obtaining an accurate approxi- 
mation of u(x,t),h(x,t) = h — h. Generally, boundary 
layer analysis deals with those problems in which pertur- 
bations are operative over very narrow regions where the 
dependent variables undergo very rapid changes across 
them. These narrow regions (sharp valley layers) fre- 
quently adjoin the boundaries of the domain of interest, 
owing the fact that a small parameter (v in the present 
problem) multiplies the highest derivative. A powerful 
method for treating the boundary layer problems is the 
method of matched asymptotic expansions. The basic 
idea underlying this method is that an approximate so- 
lution to a given problem is sought not as a single ex- 
pansion in terms of a single scale, but as two or more 
separate expansions in terms of two or more scales each 
of which is valid in part of the domain. The scales are 
chosen, so that the expansion as a whole, covers the whole 
domain of interest and the domains of validity of neigh- 
boring expansions overlap. In order to handle the rapid 
variations in the sharp valley layers, a suitable magnified 
or stretched scale and expand the functions in terms of it 
in the sharp valley regions is defined. For this purpose, 
we split u and h into a sum of inner solution near the 
sharp valleys and an outer solution away from the sharp 
valleys, and use systematic matched asymptotics to con- 
struct uniform approximation of u and h. It should be 
emphasized that at point j/o the height itself is continu- 
ous and height gradient (corresponding Burgers velocity) 
is not continuous. At these singular points the meaning 
of u± is that u±(yo, t) = u(yo±, t). Keeping in mind that 
U— > u + . the shock strength s and the shock velocity u 
are defined as s = u + — u_ and u = ^{u + + 

In appendices A and B, using the boundary layer 
method and the master equation, we have proved ana- 
lytically, that the joint moments of the height and the 
corresponding gradient difference for any a > will be 

(\6h\ a \Su\ b ) 

\x\ a+b {\Vh\ a hu\ b ) if 0< b< 1, 

\xr^({\ Vh \ a \ Vu \) + y(\s\(\u + \ a +\u-n)) if 6 = 1 



\x\"+ilp{\s\»(\u + \ a + \u-\")) 



if 1< b 



where Sh = h(xi) — h(x2), Su = u(xi) — u(x2) and 
x = x\ — X2- The quantities r\h and r\ u are the regular 
parts of d x h and d x u, respectively. For a = 0, our result 
will recover the known results for Burgers equation with 
infinitesimal viscosity [47] . 

To prove the eq.(25), we have used the fact that the 
length scale a is finite and x is let to approach zero. This 
means that we are dealing with the scaling behaviour of 
the moments (|<5/i| a |<5w| h ) for length scales \xi — X2\ << 
a. It is evident that in these length scales the height 
increments fluctuations are not intermittent. Indeed we 
find £ a = a for any moments. 



The moments of the height and height-gradient incre- 
ments i.e. (\h(xi) — h(x2)\ a ) and (|u(a;i) — u(x2)\ h ) arc 
also calculated numerically as a function of |x| = \x\ —X2\ 
for different a's and 6's. To simulate the problem . the 
KPZ equation is discretized in space and time with scales 
8 and dt respectively. The time scale dt is related to 8 
as dt = j-^-, where u m is the maximum of the height 
gradient in each time step [78,79]. At each time step the 
difference Uj = \u(xj + 1) — u(xj)\ is checked for every 
point j's. For C/j's that Uj > |<5|5, we can determine the 
positions that the height field develop a sharp valley [57]. 
Indeed this is a criterion for creation of a sharp valley in 
position Xj = yj. At points that Uj < |<5| 3 the fields 
u(xj) and h(xj) belong to the smooth part. Therefore 
the height fields h(xj) will fall into two regimes, points 
far from the sharp valleys points yj and the points in it's 
neighborhoods. For the points which the height field is 
regular or smooth, the height fields and it's correspond- 
ing gradients evolve under the KPZ and Burgers equation 
by setting the surface tension zero. Otherwise it is in the 
singularity or sharp valley region. As mentioned in the 
introduction, every sharp valley can be characterized by 
four parameters s, u, yj and h[yj). The time evolution 
of these quantities are given by the following equations 
[77] 




= au 



= -^s(h +xx - h- xx ) - f x 

d , . a , 
^ s w) = ij s ( h +xx + h- xx ) 

|%,i) = -|(4« 2 - s 2 ) + /-7 (26) 

where 7 = ht- 

To calculate numerically the scaling behaviour of mo- 
ments with x when x << a, a periodic one dimensional 
substrate consisting a discrete N-point height field with 
the length of 10000 is used. Starting with a flat initial 
condition the height and its gradient fields evolve in time. 
We consider the random force as a white in time, smooth 
and periodic in space random function which it's spatial 
correlation length is of the order of period of the given 
periodic function. To generate this type of forcing we 
use the kicking method which recently has been used in 
[51], to simulate the Burgers turbulence. The basic idea 
is that the random force can be decomposed as follows, 

f(x,t) = ^2fj(x)8(t-tj) (27) 

i 

where 8 is the Dirac distribution and where both the "im- 
pulses" fj(x) and the "kicking times" tj are prescribed 
(deterministic or random). The kicking times are or- 
dered and form a finite or infinite sequence. In this 



article the impulses are always taken smooth and act- 
ing only at length scales a. Newman & McKane [80] 
have used similar kicking, in a context where the forced 
Burgers equation is used for the study of directed poly- 
mers. Kraichnan [48] has considered a simple model in 
which there are non-smooth impulses creating directly 
saw-tooth profiles in the velocity in Burgers turbulence. 
Here the time intervals are equal to the time steps of the 
algorithm's run. 

In Figs.(l) and (2), we illustrate different snapshots 
of the time-evolution of the height, considering different 
correlation lengths a for the random periodic force, until 
the time that the system finally reaches to it's stationary 
state. The following type of kicking force is used 

F(x, t) = A(t)[cos(kx - (p{t)) + - sin(fcr - (p(t))], (28) 

3 

where A is a white Gaussian random variable in time, 
which is the noise amplitude and (p is a homogeneous 
random phase. Choosing different values for k, leads to 
different values for a. The length scale a is of the order 
of the period of F. In Fig. (3), the log- log plot of the 
moments of height increments are sketched numerically 
for a = 0.7, 5 and 7, respectively. We have found the 
exponents £ a = ra + q, where r = 1.00 ± 0.01 and q = 
—0.0012 ± 0.0002. The scaling behaviour of the moments 
of height gradients increments for length scales x « a 
is also checked. The results implies that with a good 
precision (|<5u| b ) scales with x with exponent 1 for b's 
larger than one, and scales with x with exponents £j = 
b, for b's smaller than one with precision ±0.001. The 
behaviour of £& vs b is also plotted in fig. (3). 

In summary, we study the problem of non-equilibrium 
surface growth described by the forced KPZ equation in 
1+1 dimensions. The forcing is a white in time Gaus- 
sian noise but with a Gaussian correlation in space with 
variance a. Modeling a short range correlated noise, wc 
restrict our study to the case when the correlation length 
of the forcing is much smaller than the system size. In- 
vestigating the stationary state, a general expression of 
the mixed correlations of height-difference and height- 
gradient difference at any order, in terms of the length 
scale \xi — X2 and quantities which characterize the sharp 
valley singular structures is given. Through a careful 
analysis being done over the behaviour of the sharp val- 
ley environment, wc decipher the intcrmittency exponent 
of an arbitrary a-th moment, i.e. (\h(xi) — h(x2)\ a }- 
It is proved that the height increments fluctuations are 
not intermittent and its a-th moments for length scales 
\x\ — x-2\ < <J scales as \x\ — X2^ a , where £ Q = a. In the 
present paper the limiting of v — > is taken into account 
only for finite a . Still the forcing correlation length is 
much smaller than the system size and height correlation 
length. But the limit a — > is a singular limit in our 
calculations, and moreover, it is not a priori clear that 
the limits of v — > and a — > commute at all. Us- 



ing stochastic equations which are governed over the dy- 
namics of quantities characterizing the sharp valleys we 
simulate directly the problem and check the exponents. 
We have generate the forcing using the kicking method. 
Our simulation confirm the analytical results. We believe 
that the analysis followed in this paper is quite suitable 
for the zero temperature limit in the problem of directed 
polymer in the random potential with short range corre- 
lations [81]. The same method applied to KPZ equation 
in higher dimensions would be definitely one of the con- 
sequent goals of the present work. The main message 
which might be encoded in the present work is the im- 
portance of the statistical properties of the geometrical 
singular structures for understanding the strong coupling 
regime of Kardar-Parisi-Zhang equation. 
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Appendix A 

In this appendix we are going to prove that the G-term 
in eq.(24), has a finite value in the limit v — -> 0. As shown 
in section 2 the G-term can be written as 

G = G(hi,ui,h 2 ,u 2 ,x,t) + G(h 2 , u 2 , hi, m, —x, t). (29) 

Here we prove that in the vanishing surface tension 
limit, the G term can be written as 



G(hi,ui,h 2 ,u 2 ,x,t) 



P 



ds s 



ui —s/2 



du [ui — %i)T(h\, u, s, h 2 , u 2l x, t) 



' tti+s/2 

where T(h\, u, s, h 2 , u 2 , x, t) is the PDF of 
(hi,u(yo, t), s(y ,t), h 2 (y + x, t), u 2 (y 



x,t)) 



In these intervals, boundary layer analysis can be used to 
obtain an accurate approximation of Ui(x, t) and hi(x, t). 
The basic idea is to split Ui and hi into the sum of an 
inner solution near the sharp valleys and an outer solu- 
tion away from the singular point, and using systematic 
matched asymptotic to construct uniform approximation 
of iti and hi (for details see, e.g., [77]). For the outer so- 
lution, we look for an approximation in the form of a 
series in v 



hi 

Ui 



u 



r^out 
out 



„,0 



1 



vu 



Then w° and h® satisfy 
, n ol 



(W = / 



fxi 



(32) 



i.e. Burgers and KPZ equations without the surface ten- 
sion terms. In order to deal with the inner solution 
around the singularity, let yt = yi (t) be the position of a 
shock, and define the stretched variable m = (xj — yi)jv 
and let 



l (x,t) 



y% 



s,t 



where S is a perturbation of the sharp valley position to 
be determined later. Then, <>, satisfies 



a(vi 



vi)v iz 



(33) 



where ui = dyt/dt, 7 = dS/dt and, to 0{v 2 ), vf can be 
evaluated at Xi = yi and can thus be considered as a 
function of t only. 

We study eq.(33) by regular perturbation analysis. We 
(301)3ok for a solution in the form 



Vi 



To leading order, from eq.(33) we get for v® the follow- 
ing equation 



conditional on yo being a sharp valley position. 

Let us now prove the eq.(30). Assuming spatial crgod- 
icity, for example the average of one of the terms in G, 
which is proportional to v. can be expressed as 



Ui)Vi ; 



% ZiZi 



(34) 



The boundary condition for this equation arises from the 



v(u. lXlXi \hi,h 2 ,ui 1 u 2l x)P = u(uix iXi (x,t)6(ui - Ui(xi, 
5(u 2 - u 2 (x 2 ,t))S(hi - hi(xi,t))S(h 2 - h 2 (x 2 ,t)} 

L/2 



t)) 



matching condition with u° q 



Q{v 2 



lim v\ 

Zi — >±oo 



lim 



r N 1 
v lim 

L^oo L N 



dxiUi XzXt (xi,t)8(ui -Ui(xi,t)) 



-L/2 

S(hi 



hi(x u t)). 



Clearly, in the limit as v — > only small intervals 
around the sharp valleys will contribute to the integral. 



where Sj = Sj(t) is the sharp valley strength. It is un- 
derstood that for small v matching takes place for small 
values of \x% — y%\ and large values of \z%\ = \xt — yi\jv. 
(^jDThis gives 



— tanh 



These results show that, to 0(v), eq.(31) can be esti- 
mated as 



v{uix,x i \hi,h 2 ,u 1 ,u 2 ,x)P 

v N 1 
= v urn -- 

L^oo L N 



V / dx i u i £. x .(xi,t)5(ui-u l ™(x i ,t)) 



6(hi-hr(y h t)) 

= v lim 

L^oo L N ^ 
i 

5(hi-h?{ yi ,t)) 
= v lim — — > 

i 



Appendix B 

The main aim of this appendix is to calculate the mixed 
moments (|(5/i| a |<5it| b ) by the use of master equation de- 
rived in the section 2. As we will see the term G has 
an essential role in the results being to obtain the mo- 
ments. Considering in mind equations eqs.(21) and (22), 
G 6 could be written as 

G 6 (£,u,,x,t) = 



dziu%. z .5{ui - v^{z h t)) 



dziViiziZ l m 8(iii - w") 



f 11 1 ^,11 £ , £ U W N 

/ dhduG(h — -,h + -,u — — ,u + — ,x, t) 

'-, -x,t). 



dhduG(h + — , h — —,u 



2' 



(39) 



It is proved in appendix A that the G-term can be 
(35) written as follows 



where tti is a layer centered at y\ with width 3> 0(v). 
Going to the stretched variable Zi = (xi — yi)jv, and 
using the eq.(34), we have 



G(hi,ui,h2,U2,x,t) 



dzv, 



Ozz 



dv, 



ttf) 



VOz 



advo{vo — u) 



ds 



u i — s / 2 



ui+s/2 



du (ui — u)T(hx, u, s, h>2, U2, x, t) 



»i 
(40) 



so by taking the limit as L — > oo, the z integral can be 
evaluated exactly 

v(ui XlXl \hi,h2,u 1 , u 2 , x)P 

= ap du ds T(hi, u\, Si, h 2 , u 2 , x; t) 



ill— si/2 



Ui+si/2 



dv° {v° - Ui)S(ui - v°). 



(36) 



Where (hi,ui,si,h 2 ,u 2 ,x;t) is the PDF 
of (h 1 (y 1 ,t),u 1 (y 1 ,t),s(y 1 ,t),h 2 (y 1 +x,t),u 2 (y 1 +x,t)) 
conditional on y\ being a sharp valley location and the 
spatial difference of the heights h\ and h 2 be x. Hence, 



where T(h\, u, s, h 2l u 2 , x, t) is the PDF of 

(hi, u(y , t), s(yo, t), h 2 (y + x, t), u 2 (y + x, t)) 

conditional on j/o being a sharp valley position. It should 
be emphasized that when we say y$ is a singular point, 
we mean that however the height itself is continuous at 
yo the height gradient (corresponding Burgers velocity) 
is not continues at these points. At these singular points 
the meaning of u± is that u±(yo,t) = u(±x,t)) keeping 
in mind that u_ > u+, while the singularity strength s 
and u are defined as s = u+ — u- and u = ij( u + + uJ). 
We define h + (y ,t) and h + {yo,t) as 



v{ui xixi \hi, h 2 , u 1 ,u 2 ,x)P 
r Q rm—s/2 



= —ap 



ds 



ui+s/2 



du (ui 



h+(yo,t) = h(y ) + - 
h-(yo,t) = h(yo) - | 



For late use we note that the G-term can be written in 
a more convenient manner as 

G(hi,ui,h 2 ,u 2 ,x,t) = 

f0 



(41) 

u)T(hi,u,s,h 2 ,u 2 ,x,t). (37) 

Due to the continuity of h the limit e — > is not sin- 
gular. Now let us rewrite the G s in a manner to be more 
convenient for the rest of the calculations. For this pur- 
pose let 



ds s(T(hi,ui - -,s,h 2 ,u 2 ,x,t) 



+T(h 1 ,u 1 + -,s,h 2 ,u 2 ,x,t)) 

+p / d(3 I ds sT(h\,ui + (3s, s,h 2 ,u 2 ,x,t). (38) 



Su + (x,y ,t) = u(y + \x\,t) - u+(y ,t), 
Su-(x,y Q ,t) = u-(y ,t) - u(y - \x\,t) 
5h + (x,y ,t) = h(y + \x\,t) - h + (y ,t), 
Sh-(x,y Q ,t) = h-(y ,t) - h(y a - \x\,t) 



and define U±(e, s,5h±,5u±,x,t) be the PDF's of 
(e,s(y ,t),Sh±(x,y ,t),Su±(x,y ,t)) conditional on y Q 



being a sharp valley position. Then G can be expressed 



P 5 (^LU 1 X 1 t)=8{L0)5{0+o{l). 



G\i, u, x, t) = G% (£, w, x, t) + G s _ w, t) (42) 
where 



f 

a- / dss[£/±(e, s, sgn(x)£ — e, sgn(x)u) — s, x, t) 

^ J-oo 

+U±(e, s, sgn(x)£, sgn(x)uj, x, t)] 



Define 

A(£, u, t) = lim x- 1 (P 5 (£, w, x, t) - 
= limM(Z,u>,x,t). 



(48) 



Taking the limit as x — > in the equation for P s ( 
cq.(23) ) and considering that the system has reached to 
the stationary state, it follows that A satisfies 



= -auA - 2a J dw' H(u, -u>')A(£, u>', x, t) 



ap dss d/3U±(e,s,sgn(x)£- -,sgn(x)u-/3s,x,t). (43)- B(£, u>, t) 



(49) 



Bfe,w,t) = lim G A (£,uj,x,t). 

x—*0 



ttt j . i-ii r ■ j where we have used linx r _>o(-K'(0) — K(x)) = and also 

We are interested m scaling behaviour ol mixed mo- d fi d v / v n 

ments in small length scale x. In the limit x — > it 

should be noted that P s can be decomposed into two 

parts as 

rji/s. ,s , ,, , ,, . To evaluate B note that as x — > 

.r (£, w, x, t) = Pns{x, t)P (£,u,x,t\no sharp valley) 

Su±(x,y ,t) -> 0. 

+ (I — p ns (x,t))P S (l;,uJ,x,t\sharp valley) (44) 

This implies that, as a; — > 0, 



where p ns (x,t) is the probability that there is no sharp 
valley in + x) and P 5 ^, x, £|no sharp valley) is 
the PDF of Su(x,y,t) and 5h(x,y,t) conditional on the 
property that there is no sharp valley in [y, y + x). Also 
P (£,u},x,t\sharp valley) is the PDF of 6u(x,y,t) and 
5h(x,y 7 t) conditional on the property that there is at 
least one sharp valley in [y, y + x). Since by definition of 
number density of sharp valleys p we have 



U±(8,Z,w,x,t) ^ S(8,t)6(w)8(Q 



(50) 



Pns = 1 - p\x\ + 0(x) 



(45) 



where S(s,t) is the PDF of s(yo, t) conditional on yo be- 
ing a sharp valley location. Hence, from the expression 
for G 5 , 

J9(£, u, t) = apujS(uj, t)8({) +ap<s> 5(w)5{{) 
+ 2ap5{0 dw'S{u>',t)-2apH(u)5(£) (51) 



P 5 (£, u, x, t\sharp valley) = R(£, u, x, t) + 0(1) (46) 

where R(£, s, x, t) is the PDF of £ = h(yo + x) — h(y ), 
s(yo,t) and x, conditional that yo be a shock position. 

p ns (x, t)P (£, u>, x, t\no sharp valley) 



(I- p\x\)^Q(^,- 1 t)+o(x) 

X z XX 



(47) 



here Q(r)h,rj u ,t) is the PDF of rjh(x,t) and rj u (x,t), the 
regular part of the velocity and the velocity gradient, 
respectively. Indeed we have considered the case x > 0. 
The case x < can be treated similarly. We note that, in 
the limit x — > 0, because of dealing with regular points, 
we have 

x 2 P s (xr lh ,xr/ Ul x 1 t) -> Q{ri h ,r) u ,t). 
It implies that 



where H(-) is the Heaviside function and we used 
S(s,t) = for s > since s(yo,t) < 0. Inserting this 
expression in (49), the solution of this equation is 

u, t) = (S(w) + p{s)S 1 H + pS(u, t))S(0 

Here <5 1 (^) = dS(u))/du! and to used the identity 
w5 1 (w) = — 5(uj). Using the fact that p < s >= — < r] u > 
[47], we can be restated A(£, oj, t) as 

A($,u,t) = (S{u) - (Vu)S x (cj) + pS(cJ,t))5(t). 

Hence, combining the above results, we have 

P 5 (£, u>, x, t) = (S(lj) - x(6(cj) + (nJSHu) - pS(u, t)))S(0 



o{x). 



(52) 



Which is correct for x > 0. We Reorganize this expres- 
sion as 



P s (^,u>,x,t) = [(1 — px)(8{uS) — a;(r? u )i5 1 (w)) + xpS(u>,t)]8((,) Because the singular part of Z is cancelled by g, the 

first term should be of order of o(x) 

+ o(x) 

and then we use the identity 



b rS 



dedwiei°|wiv 



5(u) -xMS 1 ^) = \q 



t +o(x). 



Now we decompose the fields h and u in terms of the 
their regular and singular parts as 

h x (x, t) = r] h (x, t) + Y^ e (yj^) S (y - Vj) 



and 



u x (x, t) = rj u (x, t) + 2J s(yj,t)6(y - yj). 



So if we let f s (iv, x, t) be defined as 

f s (^to,x,t) = (1 - p\x\)±Q (§, + ^pP^u;,*,*) 

then we can write 

P*(£,w,a;,t) = f{£,w,x,t) +o(x). (53) 

Now we can prove eq.(25) for < b < 1 and an arbi- 
trary value of a. The proof for other values of b is similar. 
Let 



(l-p\x\)-Q 
x z 

g s (£,u,x,t) = 



t) +|x|pi?(e,w,x,t), 



Because the sharp valley points have contribution in 
large a/s we can write for M > 



^iei o ki 6 (z 5 -/ < ') 



|<M 



/ ^ien W | h (z 5 -/ 5 ). 

J\u)\>M 



(54) 



oj|>M 
a+b 



M a+6 (i - Mp) 



|7J u |>M/a; 



drihdrju \rj h \ a \r) u \ h Q(r) h ,r)u,t) 



+ \x\ P / dedwlfl^wl^.w.^t) = o(x a+b ). 

Ju)>M 

We can write 

i?(£, u,x,t) as ^i?(£+, w, 1^1; *) + w : — M j t) where 

for R(£ + ,u>, \x\,t) and i?(£_, w, |x|, i), we have the condi- 
tion that £ = /i(2/o+M)-/i(z/o) and£ = h(y +\x\)-h(y ), 
respectively. When i-»0we can write £± = w±|x|. 

Since M can be made arbitrarily large, we get 

d£du\t\ a \u\ b (Z* - f) < o(x a+b ) + S m O{x) 
where 8m — * as M — > +oo. Noting that 

f Nl a+b (l^| a |^| b ) if < 6 < 1 

~ I M a+1 ((\vh\ a \vu\) + lp<M(Kl a + M a )» if b = 1 

We obtain (25) for < b < 1. For 6 > 1 the leading 
term in our calculation, will be the second term of the 
cq.(56) with the order of o(x a+1 ). The leading term is 

ip<H(Kr + M«)}. 

Also there is an alternative method to prove the eq.(25) 
for the case b > 1. The method is based on the calcu- 
lation of the mixed moment {(\h(xi) — h(x2))\ a \(u(x\) — 
u(x2))\ b ) for integer orders while b > 1, directly from the 
PDF's equation (23) by integrating over two co's, i.e. 



auP* -2a du' H(u' - w)P*(f , u', x, t) 



+2(K XX (0) - K xx {x))P^ + 2(K(0) - K(x))P^ 
+ G 5 (t,u,x,t). 



(57) 



In the limit x — > 0, keeping a finite, and in stationary 



Because of eq.(52), the first term at the rhs of eq.(54) state > [t wil1 bc simplified to 
is o(x). To estimate the second term, note that for M 
large enough 







ocjP" - 2a / du' H{u' - uj)P°{S, ,w',x,t) 



\lu\>M 



d^doj\^\ a \u\ b Z d < I d£,duj\£\ a uj 2 Z s 



u \>M 



G s (^u;,x,t) 



(58) 



< 



d^duj\i\ a uj 2 {Z s - g 



2„<5 



d£du\Z\ a uj 2 g' 



u \>M 



o(x) + \x\p 
o{x). 



d£dtu\£\ a Lj 2 S{Lj,t)8{Z) 



(55) 



First of all the term / d£du)\€\ n \w\ m G s (£,w,x,t) should 
be calculated in the x — > ±0 limit. This can be done by 
using the relation (43). Note that 



<%du;mrG s (z,u,x,t) 



a^(s\(6u + + S ga(x) S )\ m \(Sh + + sgn(x)e)\ n ) 



a^(s\Su + \ n \Sh+\ n ) 



-ap J dp (s\(Su+ + 0sgn(x)s)\ m \(5h + + sgn (x)-)\ n ) 
+ a^(s\(6u-+sgn(x)s)\ m \{5h-+sgn(x)e)\ n ) 



+ a^(s\Su-\ m \Sh-\ n ) 



-ap d(3 (s\(Su- + f3sgn(x)s)\ m \{Sh_ + sgn(a:)-)| n ). (59) 
Jo 2 

If we go back and look carefully to the definition 5u± , 
we see that Sh± ~ u±\x\ ~ o{x) and 5u± ~ o(x) as 
x — > 0. While it should be realized that the sharp valley 
strength s is of the order O(l) as e — > 0, so in the limit 
x — > 0, the result of the integral would be simplified as 

d£cM£rM m G 5 (^,x,i)~ 
a m — 1, 



2 v m + 1 



) P N" +i (| S r +i (Ki" + i^i")}. 



Finally multiplying the terms of equation (57) in |£| r 
and \uo\ m and integrating respect to £ and tu we have 



ru m+l ) 



-p\x\ n+1 (\s\ m+1 (\u+\ n +\u-\ n )) if m>neAf (60) 

where the result coincides perfectly with eq.(25) which is 
the general form of eq.(60). 
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